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Abstract 



It is known that the Standard Model does not describe all phenomena related to 
subatomic particles. This thesis presents feasibility studies of the measurement of 77 
meson decay using WASA-at-COSY detector which tests predictions of the Standard 
Model. The aim of this thesis is to estimate the time of measurement for which the 
current branching ratio upper limit of 77 — > 47r° decay can be improved. 

In order to estimate the time of measurement Monte Carlo simulations were 
performed. Kinematics of the pp — )■ ppr] — )■ ppAix^ ppSj reactions was simulated 
using PLUTO program. Next, by means of the GEANT3 program the response of 
the WASA-at-COSY detector was simulated for each particle. Received signals were 
analysed in RootSorter software package based on ROOT. 

Studies of the reaction with many gamma quanta in the exit channel required 
the investigation of WASA-at-COSY calorimeter functionality. Tests of a cluster 
building algorithm with emphasis on merging and splitting of detected signals were 
done. 

The obtained result shows that only about 3% of 77 — > 47r° decays can be properly 
reconstructed. The identification of the 77 — > 47r° — > 87 decay was based on the 
invariant mass method for both 7r° and 77 mesons. As a result, the resolution of the 
invariant mass determination of the 77 meson was estabhshed to cr = 31 MeV/c^. The 
acceptance of the WASA-at-COSY detector for pp —> ppr] — )■ ppAn'^ — )■ ppS'y reaction 
was determined to = iV^ detected/ produced = 1.5%. The time of measurement 
for obtaining a statistical precision equal to the current branching ratio upper limit 
for Theam = 1.3 GeV is about 13 hours. It has also been established that the lower 
the excess energy in the 7] production is the shorter the time of the measurement 
needed. 
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1 Introduction 



Since last century the concept of symmetry has become very important in physics. 
The 1918 year, in particular and famous Noether's theorem, which relates symme- 
tries to conservations laws is arguably considered as a crucial moment in the history 
of the concept of symmetry. The first one who regarded symmetry principle as the 
primary feature of nature was Albert Einstein. He recognized the symmetry implic- 
itly in Maxwell's equations and elevated it to a symmetry of space-time itself [1]. 

The symmetries have found applications in the Standard Model (SM) which 
is a theory describing the electromagnetic, weak, and strong interactions, which 
govern the dynamics of the known subatomic particles. However, although almost all 
experimental data can be explained by the SM, it does not explain basic phenomena 
such as the existence of three families of fundamental fermions, the origin of CP 
violation, etc. It is expected that the SM breaks down at some point. A place where 
New Physics can be looked for is the limit of validity of the basic symmetries of 
charge conjugation (C), parity (P) and time reversal (T), as well as combination of 
these symmetries CP and CPT 0. 

In the Standard Model CP violation is described by the phase in the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix. Six quark flavours are grouped into three 
families. CP violation is related to family- changing interactions, while in family- 
conserving cases CP violation is not included in the SM. It is deemed that detailed 
studies of CP violation may lead us to New Physics that goes beyond the Standard 
Model. Decays which test flavour-conserving P as well as CP violation are t] and rj' 
decays into 2n [2J. 

1] meson is a part of the light pseudoscalar nonet together with mesons vr, K and 
f]'. Real f] and rj' mesons are a combination of the octet and singlet states {r]s and 
rji) of the SU(3) symmetry. The t] meson is eigenstate of C, P and CP operators, 
Cirf) = +1, Pirf) = —1. Therefore, the decays of the rj meson constitutes a tool to 
study these symmetries. The decay of r/ — )■ 27r° is the simplest example of a decay 
violating P and CP symmetry. The branching ratio for this decay is predicted to be 
less than 6 -10"^^ [3] and a present experimental upper limit amounts to 3.5-10"^ [1]. 

CP symmetry test 1] 27r° is hard to improve due to the background connected 
with direct 27r° production in every t] production reaction. A possible new test is 
the decay into four pions [2]: 

1] 47r° 

The present upper limit for branching ratio BRij] — )■ Atx^) < 6.9 ■ 10^'' was 
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determined by Crystal Ball Collaboration [5]. 

CP symmetry in r/ — )■ 4tt^ can be conserved if wave function of pions in final 
state involves higher partial waves. However, in this case, estimated branching ratio 
of the CP conserving decay is less than 10~^° [3]. A measurement of a non-zero BR, 
for this decay, above estimated limit would be a signal for a CP violating process. 

The main goal of this thesis is to estimate time of measurement, with WASA-at-COSY 
detector, for the reaction pp — )■ ppr] — )■ ppAn^ — )■ ppSj needed to obtain precision 
better than a present upper limit of BR{ri — > Atx^). 

Next section of this thesis comprises short theory introduction of symmetries in 
general with emphasis on parity and charge conjugation symmetry. The CP sym- 
metry violation in rj decays [r] — )■ 27r° and r] — > 47r°) are described in section [3l 
In section H] the cross section for the r] meson production as well as the multip- 
ion production in proton-proton collision are presented. Further on electromagnetic 
calorimeter and other WASA-at-COSY detectors are described in section O The 
section presents geometrical acceptance of the calorimeter and its efficiency for 
gamma quanta reconstruction. Level of statistical uncertainty of the branching ra- 
tio, achievable with WASA-at-COSY detector is determined as a function of time in 
section [71 Finally, the time of measurement needed to achieve a required accuracy 
in the studies of r/ — )■ Att^ decay is estimated as a function of excess energy. 

The work is supplemented with Appendix where results of kinematical fit used 
for pp — )• ppr] — )■ pp4:TT^ — )■ pp87 reaction are presented. 



2 Symmetries 



Symmetries can be divided into continuous and discrete: 

• Continuous symmetries lead to conservation laws: 

1. invariance under spatial transformation leads to momentum conservation 

2. invariance under time transformation leads to energy conservation 

3. invariance under spatial rotation leads to angular momentum conserva- 
tion 

4. invariance under rotation in isospin space leads to isospin conservation 

5. invariance under a change in phase of a wave function of a charged particle 
leads to electric charge conservation 

• Discrete symmetries such as e.g. 

1. Parity symmetry 

2. Charge conjugation symmetry 

3. Time reversal symmetry 

do not lead to new conserved quantities in classical mechanics. They are im- 
portant in quantum mechanics. 

Symmetries tell us which interaction is responsible for particular process. Elec- 
tromagnetic and strong interactions conserve parity and charge conjugation symme- 
try. However, it was observed that weak interactions violate these symmetries. This 
violation was first discovered in 1957 for parity symmetry [6] and it was thought 
that a combination of C and P (a CP symmetry) is conserved. Yet, in 1964 it was 
discovered in experiment with kaons, that CP symmetry is also violated. In the 
decays of neutral kaons it occurs at a level of 10~^ [7j. Next assumption was that 
a combination of C, P and T (a CPT symmetry) is conserved. Until now nobody 
has discovered violation of this symmetry. In this thesis I concentrate on the CP 
symmetry therefore P and C symmetries will be described more detailed. 

2.1 Parity symmetry 

The contents of this and following subsection are mostly based on Ref. [S]. 
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Symmetries 



Parity is related to spatial reflection. Spatial reflection of polar vector x gives 
—X. Thus the action of the parity operator P on a single-particle state represented 
by wave function ■0(r, t) gives: 

P'ilj(x,t) ^ X ■'ijj(-x,t), (2.1) 

where X is a phase factor. If wave function ■0(r, t) is eigenstate of parity operator 
then: 

Pi;{x,t)^P-ij{x,t), (2.2) 

where the eigenvalue P = ±1 is called the intrinsic parity or parity of the state. 
Parity is a multiplicative quantum number. 

If a particle has an orbital angular momentum I, there is an additional contribu- 
tion to the total parity of the system. In this case particle's wave function consists 
of a radial part Rni{r) and an angular part Yi^{9, 0): 

i^nlm{x)^Rnl{f^-Yr{e,(t>), (2.3) 

where n is principal and m is magnetic quantum number. Yi^{9,(f)) is a spherical 
harmonic. It can be shown that the parity transformation x — > —x for spherical 
coordinates implies that: 

r ^r, ^n-e, (f) ^ n + (f), (2.4) 

and 

vne, <t>) ^ rn-K - ^, TT + 0) = (-1)' . Yr{e, <t>). (2.5) 

Therefore, 

Pi^nlmix) = P ■ ^nlm{-x) = P ' ("1)' " i^nlm{x) , (2.6) 

where ipnim{x) is an eigenstate of parity operator P with eigenvalue P ■ (—1)'. 

Strong and electromagnetic interactions are invariant under parity. Parity trans- 
formation does not change the Hamiltonian of the system governed by strong and 
electromagnetic interactions. 

2.2 Charge conjugation symmetry 

Charge conjugation is the operation that changes particles into their antipar- 
ticles. The result of action of charge conjugation operator for particles which are 
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identical to their antiparticles (a = a) is different from the result if particle differs 
from its antiparticle {b ^b): 

C\a,^a) = Ca\a,iJa) (2.7) 

C\b,i:,) = \b,ij-,), (2.8) 

where = ±1 is a phase factor, letter a in ket denotes particle and tpa corresponds 
to its wave function. Phase factor for neutral pion is equal to +1 (Ca(vr°) = +1). 

Particles can be divided into two groups. The first group comprises particles 
which wave functions are eigenfunctions of C operator (12.71) . These particles are also 
theirs own antiparticles. Neutral pions, t] mesons and gamma quanta are examples of 
this kind of particles. In the second group particles and corresponding antiparticles 
are different (12. 8p . Examples for this case are charged pions and protons. 

Charge conjugation is multiplicative quantum number [9j, operator C acting on 
a system of two neutral pions ip = ipai^a gives: 

C I'ip) = Ca |a, ipa) Ca \a, ipa) = \a, ipa) \a, ipa) ■ (2.9) 
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3 CP symmetry violation in the 77 meson decay 



To understand why C and P symmetries are violated in the rj meson decays, it is 
necessary to know a wave function of the final state system. It is worth to remember 
that the total wave function of a final state is a product of spin, isospin and spatial 
parts: 



Y total Yspin ' Y isospin ' Y spatial- \^'^) 

Therefore, resultant sign of the wave function under the parity transformation is the 
product of parities of those three components. 

In the following subsections each part of the total wave function will be consid- 
ered for two and four pion systems. 



3.1 System of two neutral pions (77 — tt^tt^) 

At the beginning of this subsection spin part will be considered. Spin of neutral 
pion is equal to zero, 5'(7r°) = 0. Therefore, the spin of a pair of pions is also equal 
to zero, ^(Tr^Tr^) = 0. In the spin space the parity operation is equivalent to the 
exchange of the spins in the wave function. Using the property of Clebsh-Gordan 
coefficients [TU]: 

< simiS2m2\SM >= (-1)'^+'^-^ < s2m2Simi\SM >, (3.2) 

where Si denotes spin of i-th particle, is spin projection, 5* stands for spin of 2tt 
pair and M corresponds to its projection, it can be shown that the sign of the spin 
part under the parity transformation of the wave function is positive {ipspin = +!)• 
It is because (-l)0+o-o = +i. 

The isospin of a tt meson is equal to 1, /(vr) = 1. Hence, the isospin of two pions 
can be equal to 12-,^ = 0, 1, 2. However, the Clebsh-Gordan coefficient for two neutral 
pions in the final state of isospin 1 equals to zero. That is why the isospin of (vr^vr'^) 
system can be only I2n0 = 0, 2. Using equation (13. 2p for isospin, there are only two 
possibilities for (— l)*i+*2--'^ coefficient: (—1)^+^^° or (—1)^+^"^; in both cases this 
factor is equal to +1. 

Finally, as regards spatial part of the wave function based on equation (12.61) it 
can be written that the parity of a wave function for two particles is the product of 
internal parities of each particle and the parity resulting from their orbital angular 
momentum L: 
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CP symmetry violation in the rj meson decay 



PV(= = P(^i)P(^2)(-l)^ (3.3) 

The internal parity of pion is P(7r) = —1. From equation (13.31) one can see that 
the change of sign of spatial part of 271^ wave function under the parity transforma- 
tion is equal to the factor (—1)^. 

In order to evaluate value of L, information about the total angular momentum 
J of the decay t] — tt^tt^ can be used. Angular momentum is described by equation: 

J = L + S, (3.4) 

where L denotes orbital angular momentum and S stands for the total spin of the 
system. Because a total angular momentum of initial and final state has to be equal 
and J{ri) = so J(7r°7r°) also has to be equal to zero. Knowing that ^(Tr^Tr^) = and 
using equation (13. 4p . one can get = L+0. Therefore, the orbital angular momentum 
of n^TT^ system originating from the f] meson decay equals to zero, L(7r°7r°) = 0. Thus, 
the factor (—1)^ in equation (13. 3p is equal to +1. 

Parity symmetry in 77 — )■ n^n^ decay can be determined by gathering above 
results. From the above considerations we may infer that the 7r°7r° system produced 
in ?7 — )• 7r°7r° decay possesses parity equal to +1 because the change of sign of spin, 
isospin and spatial parts of wave function under the parity transformation are all 
equal to +1 (P(V'spm) = +1, P{i^isospin) = +1, Pii^spatiai) = +1). Taking into 
account that P(?7) = —1, one can see that P symmetry is violated in r/ — )■ 7r°7r° 
decay. 

Charge conjugation symmetry phase factor Ca of 2tt^ system with L = has the 
same value as for rj and is equal to +1. Thus, C symmetry in r/ — t- ti^tt^ decay is 
conserved, but CP symmetry is violated because P is violated. 

Consequently, in the rj — t- tt^tt^ decay a C symmetry is conserved but P and CP 
symmetries are violated. 

3.2 System of four neutral pions (r/ — > 47r'^) 

Equation (12.91) implies that charge conjugation eigenvalue of the 47r° system is 
just equal to (+1)^ = +1. 

Spin and isospin parts of total wave function are analogous to those in the 2tt^ 
system and they do not change its sign under the parity transformation. Spatial 
part of four particles system is more complex than for two and will be considered in 
detail in the following part. 

System of Air^ can be described, for example, in a way shown in Fig. 13. 1[ In this 
case there are three orbital angular momenta li, I2, L, which need to be considered. 
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Figure 3.1: Picture illustrates angular momenta in four neutral pions system. 



Due to Bose symmetry li and I2 are even [HI [ID] • Thus, in general, there are three 
possibilities of orbital angular momenta in the system: 

1. Orbital angular momenta in two pion pairs are equal, li = l2- Four neutral 
pions system is a bosonic system, so due to Bose symmetry L is also even. 
Substituting it in (12.61) . when all orbital angular momenta are even, one can 
calculate that parity of the system is equal to P = +1. Thus, in this case 
CP(47rO) = +1 while CP{r]) = -1, and so the CP is violated 

2. Orbital angular momentum of one pion pair is zero, Zi = and of the other 
is even, I2 = 2n, where n is a natural number. In this case L is even because 
J{ri) = 0, so 7(47?") also has to be zero. It is possible if and only if L and I2 
are equal. In this case CP is also violated. 

3. Orbital angular momentum of one pion pair is at least two, /i > 2 and of the 
other is at least four, I2 > 4. Considering li = 2 and I2 = 4, these two vectors 
can be combined and one of possible result is three, li^2 = 3. If L is also equal 
to three (L = 3) then the total angular momentum in the 77 — )■ Att^ may be 
conserved, however one obtains P = — 1. It is because (—1)2+4+3 _ ^^lis 
case CP{ri) = CP(47r^) = —1 and CP symmetry is conserved. 

To check how probable is this case estimation of orbital angular momentum 
values accessible in the system is needed. Excess of mass, which can be trans- 
formed into energy in rj ^ Att^ decay is approximately 8 MeV/c^. Angular 
momentum in classical mechanics is given hj L = f x p. In order to estimate 
a maximum angular momentum accessible in a Att^ system, one can make a 
conservative assumption that one pion pair do not take any energy. Assum- 
ing that remaining two pions share accessible energy equally, one can find 
|p| = 68 MeV/ c. For r it can be assumed that radius of r] meson is not larger 
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than radius of nucleoli (~ 1 fm). Substituting these values into equation for 
angular momentum we obtain L = 68 MeV ■ fm/c. The obtained value of L is 
much smaller compared to the unit of angular momentum h ~ 197MeV ■ fm/c 
Since we need at least four units in one pion pair, L is only 9% of needed 
value. In addition, it is worth to remember that to the CP symmetry con- 
servation also two units of orbital angular momentum in the second pair of 
pions are necessary. It is estimated that branching ratio for process with CP 
symmetry conservation in 77 — )• Air^ decay is less than 10"^'' [3]. 



4 Reaction pp — > ppr] and pp — > pp multipion 



The Cooler Synchrotron in Forschungszentrum Jiihch (see subsection 15. 11) enables 
rj production in two main reactions: pp — >■ pprj and dp — >■ ^He 7]. Cross sections as a 
function of ^/s for these reactions are shown in Fig. 14.11 
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Figure 4.1: Cross sections for the t] meson production in proton-proton (left) and 
proton-deuteron collisions (right). The left hand side histogram contains experimen- 
tal data from Refs. [TT], [121 [131 HH [El [IS] and the right hand side histogram shows 
data from Ref. [IHIIH]. 



One can see that the shape of the cross sections as a function of ^/s is quite 
similar. Starting close to threshold of the reaction they can be described by rapid 
increase in cross section value and later transfers into plateau. The plateau for 
proton-proton collision is at a level of 10^ fib and for deuteron-proton collision it is 
located at above 0.4 fib. Since the t] — )■ 47r° is a rare decay, a high statistics of r] 
production is required. Thus, proton-proton collision has been chosen and will be 
considered in further analysis. 

For the estimation of the precision for the extraction of a BR of 77 — )■ 47r° decays 
we need to estimate also the physical background. The background is related to 
techniques which are used for decay particles reconstruction. 
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Reaction pp — )• ppr] and pp — )• pp multipion 



At the WASA-at-COSY experiment the missing mass and invariant mass tech- 
niques are used for rj meson identification via the pp — > ppX reaction where X 
denotes the unobserved particle. If protons four-momenta P = {E,p) before and 
after the colhsion can be determined, then conservation of a total four-momentum 
vector gives: 

neam + ^target = Pi + Ps + ^X, (4.1) 

where Fx stands for the four-momentum of unobserved particle. Employing equation 
(14.1 p the mass rrix can be calculated as: 

ml = |PxP = l^eam + ^target ~ ^1 - FsP- (4.2) 

Another useful method of unstable particle identification is invariant mass tech- 
nique. It is based on a reconstruction of a rest mass of decaying particle which can be 
calculated if information about all decayed products is known. The mass is described 
by: 

y i i 

where Ei and Pi denote energy and momentum of i-th particle respectively. 

In spite of using both techniques the rj — )■ An^ decay can be misidentified due to 
the unavoidable physical background as well as due to the misinterpretation of the 
detector signals. The main physical background of — > 47r° is prompt 47r° production 
in proton-proton collision. This reaction is difficult to reject because all particles in 
exit channel for both reactions are the same. Difference between mass of rj and 47r° is 
only approximately 8 MeV/c^. Resolution of WASA-at-COSY detector is not good 
enough to reduce the background events using invariant mass technique. 

The current available data of cross section for multipion production in proton- 
proton collision are shown in Fig. 14. 2[ The cross section in the figure is shown as a 
function of excess energy with respect to the t] meson production threshold. 

The cross section for An^ production in proton-proton collisions is not established 
and so far only an upper limit for a single energy point at ^/s = 2.498 GeV/c^ was 
determined to 5 nb |22]. In order to estimate the background from direct pions 
production {pp — > ppAn^) for the pp — > ppr] — > ppAir^ reaction we have derived an 
energy dependence of the upper limit of the total cross section under assumption of 
the homogeneous phase-space population [TU] : 



a(g) = const ■(Q/v^)(3™-5)/2, 



(4.4) 
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Figure 4.2: Experimental cross section for direct pions production in proton-proton 
collision (points). Superimposed lines denotes result of the fit of formula l44l to the 
experimental points. 

where Q and s denote the excess energy and the total energy in the centre-of-mass 
system respectively, m stands for the number of particles in the exit channel. 

The result is shown in Fig. 14.21 as the lowest, green line for which the absolute 
value was fitted to the experimental upper limit. For the comparison in Fig. 14.21 
also cross sections for tt , 27r and Svr production are shown. The superimposed lines 
indicate result of the fit of formula 14.41 with amplitude as the only free parameter. 
The obtained values of const are listed in Tab. 14.11 



Number of neutral pions 


Normalization const from equation |4.4| 


1 


10^ 


2 


3- 10^ 


3 


3-105 


4 


2.35-10^ 



Table 4.1: Values of const parameter determined by fitting parameterization 14.41 to 
experimental data shown in Fig. 14.21 



In Fig. 14.21 one can see that e.g. at 50 MeV above the t] meson production 
threshold the cross section for 47r production is by about five orders of magnitude 
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smaller than for 27r production. Consequently, the physical background for studies 
of 77 — > 27r° decay is by five orders of magnitude larger than for 77 — > 47r° decay. 
Moreover, with decreasing excess energy this ratio is increasing. 



5 Experimental facility 



5.1 Cooler Synchrotron COSY 

The contents of this and the following subsection are mostly based on Refs. [231 

The COoler SYnchrotron (COSY) is located at the Institute of Nuclear Physics 
of the Research Centre Jiilich in Germany. It is a storage ring of 184 m circumference. 
The facility provides unpolarized and polarized proton and deuteron beams in the 
momentum range from 0.3 GeV/c up to 3.7 GeV/c. To decrease the momentum and 
spatial spread of the COSY beam, electron and stochastic cooling is used. 

5.2 WAS A detector 



The Wide Angle Shower Apparatus (WASA) is approximately An angle detector 
and is located at the COSY ring. It consists of two main parts: The Forward Detector 




Central Detector Forward Detector 



Figure 5.1: Schematic view of the WASA-at-COSY detector. The abbreviations of 
the name of detectors are explained in the text. 

5.2.1 The Forward Detector 

The Forward Detector measures charged target-recoil particles and scattered 
projectiles. It covers scattering angles between 3° and 18°, and consists of 30 layers 
of detectors: 



21 



22 



Experimental facility 



The Forward Window Counters (FWC) 

The FWC is the first detector in the FD along the beam direction. It consists of 
two layers and 48 elements altogether. Each element is made of plastic scintilator 
with thickness of 3 mm. The FWC is used as a trigger. Coincident hits in different 
subdetectors at the same azimuthal angle are used to select events originating from 
the interaction vertex. 

The Forward Proportional Chambers (FPC) 
The FPC is mounted directly behind the FWC tracking device downstream the 
beam. It consists of four modules, each with four layers of 122 proportional drift 
tubes (straws). The modules are rotated by 45° with respect to each other (in the 
plane perpendicular to the beam axis). The FPC is used to reconstruct tracks of 
charged particles. 

The Forward Trigger Hodoscope (FTH) 
The FTH is third sub-detector. It is made of three layers built out of 5 mm plas- 
tic scintilators. The first layer consists of 48 radial elements with shape similar to 
triangles. Second and third layers are divided into 24 elements shaped as Archi- 
median spirals. It provides hits multiplicity and both angles: polar and azimuthal. 
The FTH is mainly used in the first level trigger logic in a coincidence with other 
sub-detectors. 

The Forward Range Hodoscope (FRH) 
The FRH consists of 5 layers made of plastic scintilator, each layer is comprised of 24 
elements. First three layers have thickness of 110 mm, while thickness of fourth and 
fifth layer is equal to 150 mm. The FRH is used for kinetic energy reconstruction 
and particle identification. 

The Forward Range Intermediate Hodoscope (FRI) 
The FRI is two layered detector placed between the second and the third layer of 
the Forward Range Hodoscope. Each of its layer is made of 32 rectangular elements. 
Bars in one layer are oriented horizontally and vertically in the other. The purpose 
of placing FRI was to improve two-dimensional position sensitivity within the FRH. 

The Forward Veto Hodoscope (FVH) 
The last active element of the Forward Detector is FVH. It is made of 12 horizontal 
plastic scintilator bars with thickness of 2 cm. It is used for detecting high energetic 
particles. Lately second layer of FVH was installed. It will be used as a stop detector 
for Time-of-Flight measurements. This will improve the energy resolution for high 
energy particles [27j . 

The Forward Range Absorber (FRA) 
The FRA is a passive absorber layer made from iron plates with thickness of 5 mm. 
Maximum thickness of the absorber is 100 mm. It can be installed between the 



23 



FRH and FVH detectors. For adequate thickness of the absorber slower protons 
from reactions in which rj is produced are stopped in the absorber, while faster 
protons and pions [r] is not produced) can give signals in FVH. 

5.2.2 The Central Detector 

The Central Detector (CD) surrounds the interaction point. It consist of the 
Mini Drift Chamber, the Plastic Scintillator Barrel, the Superconducting Solenoid, 
the Scintillating Electromagnetic Calorimeter. All these components allow for mea- 
surement of momentum and energy of neutral and charged particles. 

The Mini Drift Chamber (MDC) 
The MDC is placed around the beam pipe. It is assembled from 17 cylindrical layers 
with 1738 straw tubes. It covers scattering angles between 24° and 159°. The MDC 
is used for the determination of the charged particle momenta and the interaction 
vertex. 

The Plastic Scintillator Barrel (PSB) 

The PSB is a cylindrical detector surrounding the drift chamber. It consists of a 
central barrel part (cylindrical) and two end caps (forward and backward). Caps are 
made of trapezoidally shaped 8 mm thick plastic scintillator elements. Cylindrical 
part consists of 2 layers composed of rectangular strips with a small overlap between 
neighbouring elements. The PSB is used to distinguish charged form neutral particles 
and serves as AE detector for charged particles identification. 

The Superconducting Solenoid (SCS) 
Magnetic field is necessary to determine momenta of charged particles with the 
MDC. It can be done by dint of the SCS. The solenoid surrounds the Mini Drift 
Chamber and the Plastic Scintillator Barrel. It provides axial magnetic field used for 
bending charged particle trajectories. Cooling of the solenoid is provided by liquid 
Helium. The maximum axial magnetic field is up to 1.3 T in the interaction region. 

The Scintillating Electromagnetic Calorimeter (SEC) 
The SEC enables measurement of scattering angles as well as energy of gamma 
quanta, electrons and positrons. Energy of gamma quantum can be measured in the 
range from 2 MeV up to 800 MeV. The calorimeter consists of 1012 sodium-doped 
Csl scintillating crystals grouped into three parts shown in Fig. 15.21 with different 
colours. Starting from left: forward, central and backward part. Calorimeter covers 
full azimuthal angle (0) and polar angle (6) from 20° to 169° which corresponds to 
96% of full solid angle. The crystals are shaped as truncated pyramids to conform 
spherical geometry of the calorimeter. Three parts of the calorimeter differ in polar 
angle range and crystal size. The lengths of the crystals vary from 30 cm (16.2 
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Experimental facility 




Figure 5.2: The Scintillating Electromagnetic Calorimeter (SEC) used in detector 
WASA-at-COSY. It is divided into three parts, marked in different colours, starting 
from left: forward, central and backward part. 

radiation lengths) to 20 cm. The longest are in the forward part, the shortest in the 
backward part. The forward part consists of 4 layers with 34 elements each covering 
polar angle from 20° to 36°. In the central part 9 is between 34° and 150°. It consists 
of 17 layers and 48 elements in each layer. The backward part covers polar angle 
from 150° to 169° and is built out of 3 layers. Two of them has 24 elements and one 
layer possesses only 12 elements. 

Size of crystals depend on part in which they are placed and are crucial for 
angular resolution of the calorimeter (see subsection 16.31) . The energy resolution for 
gamma quanta is described as ^ = y - -^^^ stopped charged particles the 

energy resolution is equal to about 3%. 

Each crystal is read out individually by photomultipher. Readout is done outside 
of the iron yoke in order to avoid influence of magnetic field. 

More details on the construction and the design of the calorimeter can be found 
in Ref. |25|. 



6 Identification of 77 — An decay 



6.1 Geometrical acceptance of calorimeter 

The first step of reaction pp — )■ pprj — )■ ppAn^ ppS'y analysis is to check what 
fraction of all reactions can be fully observed with the WASA detector. For this 
purpose 1 000 000 reactions of r] meson production and decay into four neutral 
pions with the homogeneous phase-space distribution were generated by means of 
PLUTO program [28]. Output file from the PLUTO was used in the WASA Monte 
Carlo (WMC) detector simulation. The WMC is based on the GEANT3 (GEom- 
etry ANd Tracking) software package developed at CERN [21]. The GEANT pro- 
gram describes the passage of particles through the matter and simulates response 
of WASA-at-COSY detector. The analysis of detector signals is carried out using 
RootSorter. The RootSorter is a software package based on ROOT [30] . 
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Figure 6.1: The number of reconstructed clusters for pp — )■ ppr] — )■ ppAir^ — )■ ppS'y 
reactions as a function of the number of gamma quanta which escaped from the 
Scintillating Electromagnetic Calorimeter. Histogram is drawn under the condition 
that two charged particles are reconstructed in the Forward Detector and there is 
no charged particle in the Central Detector. 

The result of that studies is shown in Fig. 16. 1[ The figure presents the number of 
reconstructed clusters as a function of the number of gamma quanta which escaped 
from the Electromagnetic Calorimeter. A gamma quantum can escape if it flies into 
the beam pipe or the pellet pipe. Histogram is drawn under the condition that two 
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charged particles were reconstructed in the Forward Detector and no charged particle 
in the Central Detector. These conditions are imposed in order to choose reactions 
when two protons fly into the FD and to reduce events when electron-positron pairs 
are created in the CD material. Reconstructed signal in SEC is considered to be 
a charged particle if corresponding signal in PSB can be found. Only 38.7% of all 
simulated events obey the above conditions. 

When one or more gamma quanta escape from the SEC the information which 
they carry is lost. Therefore, for further consideration, only events when all eight 
gamma quanta hit the SEC are taken into account. Distribution of reconstructed 
clusters for the situation when all eight gamma quanta hit the SEC is shown in Fig. 
16.21 It reduces number of accepted events to 11.2%. 
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Figure 6.2: Projection of Fig. 16. li on X-axis for events when all eight gamma quanta 
from reaction pp — )■ ppr] — )■ ppAir^ — )■ ppSj hit the Scintillating Electromagnetic 
Calorimeter. 

The histogram 16.21 is plotted for events where no gamma quanta escaped from 
the SEC. One can see that the number of identified cluster varies from 3 to 10. This 
is a manifestation of merging and splitting of detected signals which leads to wrong 
clusters reconstruction (see subsection 16. 3p . 



6.2 Pions reconstruction 

For reaction where eight gamma quanta are in exit channel, in the further step 
of the analysis it must be established which gamma quanta originate from a decay 
of the same pion. That is why a routine which matches each two gamma quanta into 
pairs was elaborated. Purpose of the routine is to identify pairs of gamma quanta 
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originating from the decay of the same pion. The prepared procedure first calculates 
squared invariant mass of every possible pair of gamma quanta and subtracts the 
squared mass of 7r°. Next the results are combined into sets of four pairs which in- 
cludes all eight gamma quanta. Further on absolute values of the differences between 
the invariant mass of the pair and vr" mass within a given set are added together. The 
set for which the sum of differences from the pion mass is the smallest is assumed 
as a correct assignment and is taken for the further analysis. 

If a number of reconstructed clusters is less than eight, there is not enough 
information to reconstruct the reaction pp — )■ pprj — )■ ppAir^ — )• ^^87. That is why 
reconstruction of the rj meson decay is done only for events when eight or more 
clusters are observed which reduces the number of accepted events further from 
11.2% to 4.8%. 




Figure 6.3: Invariant mass of pairs of gamma quanta originating from the rj meson 
decay, selected by means of algorithm described in the text. The left-side histogram 
corresponds to events with eight reconstructed clusters, the right-side histogram 
contains events with eight or more reconstructed clusters. Gaussian functions (16. ip . 
fitted to the histograms, are shown as solid line. 

Fig. 16.31 shows invariant mass of pairs of gamma quanta identified by above 
described algorithm. Left panel presents the result for events with eight clusters 
reconstructed and the right histogram shows events when eight and more clusters 
have been reconstructed. To both histograms Gaussian function is fitted: 



28 



Identification of r/ — )• Att^ decay 



^ f — ix — Meanr\ 

Constant ■ expl — . (6.1) 

V 2 ■ bigma^ / 

Values of fit parameters are included inside the figures. Mean and sigma values do 
not differ significantly. However, statistics is larger for events with eight and more 
clusters reconstructed. Therefore, these events will be considered during further 
analysis. 

Invariant mass spectra in Fig. 16.31 are not symmetric and one can see tails at the 
left sides of both histograms. The tail is related to the functioning of two procedures. 
The first one matches gamma quanta into pairs and is described above. The second 
is algorithm which bases on energy deposited in the SEC and changes it into clusters 
information. To understand this process, knowledge of the cluster building routine 
is needed. 

6.3 Cluster building routine 

Description of the cluster building routine is mostly based on Refs. [2S1 EI] • 
Gamma quanta give signals in the SEC crystals by electromagnetic showers pro- 
duction. Charge of the signal provides information about energy deposited in a given 
detection module. Every crystal has its own readout system so it is taken as a sep- 
arate element. To reconstruct observed particles, the routine combines all available 
hits into clusters. 




Figure 6.4: Ilustration of cluster creation. 

Cluster finding routine is an iterative procedure. First, it searches for the module 
with the highest energy deposition. The module becomes a reference element for a 
new cluster. Later it creates a square of 3 x 3 elements with the reference element 
in the centre. If any element with deposited energy is inside the square, it will 
be added to the cluster. All new added elements become also reference elements 
in next step of procedure. The procedure stops reconstruction of a cluster when 
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no new elements with an energy deposition higher than 2 MeV in the square are 
found. Further on the routine repeats the cluster finding excluding elements assigned 
to the already identified clusters. The procedure is repeated until all modules with 
deposited energy are attached to clusters. An example of cluster with structure of hit 
crystals is shown in Fig. 16.41 In the shown example elements with lined background 
are parts of a reconstructed cluster. One element, at right side is not attached to 
the cluster. 

Restrictions for neighbouring elements which can be added in a cluster are in- 
cluded in the routine. Default maximal time difference between elements must be 
smaller than 50 ns. Minimal energy deposits are 5 MeV for a new cluster starting 
element and 2 MeV for a neighbouring element. Total energy deposited in a cluster 
is a sum of energy deposition in all crystals from which the cluster is built. Total 
energy of a cluster must be larger than 10 MeV. If it is smaller than 10 MeV, the 
cluster is rejected from further analysis. 

The cluster time and cluster coordinates are calculated as a mean value of all 
attached elements weighted by the energy deposition. Energy deposited in a cluster is 
used to calculate gamma quantum kinetic energy. Particle identification as a neutral 
requires no corresponding signal in the PSB. 
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Figure 6.5: Probability for reconstruction of a given number of clusters for two 
gamma quanta hitting the calorimeter as a function of opening angle between them. 
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Test of the ability of the used algorithm to reconstruct clusters depending on the 
opening angle between two gamma quanta was performed. Two gamma quanta with 
different opening angle were simulated in the test. Results are shown in Fig. 16.51 If 
opening angle between two gamma quanta is small, they are recognize as one cluster. 
It is enough that there is no 'empty' crystals between elements with non-zero energy 
deposition information. One cluster reconstruction dominate for small value of the 
opening angle (red line). Value of opening angle for which two cluster recognition 
is starting to be dominant depends on part of the calorimeter. It is connected with 
different size of crystals in different part of the SEC. 
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Figure 6.6: Histogram of energy deposited in an additional cluster (splitoff) as a 
function of opening angle between cluster generated by gamma quantum and the 
additional cluster. 

In Fig. 16.51 apart from merged signals reconstructed as a single cluster (red line) 
one can see events when three clusters were recognized (blue line). Reconstruction 
of three clusters results from splitting of signals. One of possible reason for splitting 
is that electromagnetic shower caused by gamma quantum may spread across few 
modules. The middle crystal can be omitted in case if deposited energy is less than 
2 MeV or gamma quantum do not react in it at all. Thus, module read by the 
routine as 'empty' will be produced. To check this hypothesis, energy deposited in the 
additional cluster (splitoff) as a function of opening angle between cluster generated 
by gamma quantum and additional cluster is shown in Fig. 16.61 No entries below 
10 MeV are due to the requirement that the total energy of every cluster must be 
larger than 10 MeV. Distribution in Fig. 16.61 shows that splitoffs are characterized 
by small energy deposition and small distance to real particle. But there is a small 
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fraction of events for which the spht cluster is far from the hit position of the gamma 
quantum. 

6.4 T] meson reconstruction 

At this point, when we aheady know how cluster building routine is working, we 
can go back to explain tails in the histograms shown in Fig. 16.31 As it is mentioned 
above, tails in histograms (16.31) are artefacts of used algorithm and the granularity 
of the calorimeter. 

It may happen that in a given event number of merged signals is equal to the 
number of split ones. In this case one will reconstruct eight clusters as expected for 
the pp — )■ ppt] — )■ ppAn^ — )• ppS'y reaction. Such event will fulfil all required criteria 
but nevertheless it will lead to the wrong reconstruction of energy and momenta of 
some gamma quanta. It manifests itself as a tail in the invariant mass histogram. 
In order to demonstrate this effect in Fig. 16.71 invariant mass for pairs composed of 
real cluster and the additional one are plotted. 




Figure 6.7: Histogram of invariant mass of cluster generated by gamma quantum 
and an additional one (splitoff). 

To held as little events with splitoffs as possible, only these pairs whose invariant 
mass differs from the pion mass by no more than 2a are used in the further analysis. 

Events in which all four gamma quantum pairs fulfil restrictions mentioned above 
are used for the t] meson reconstruction. Invariant mass of all eight gamma quanta 
is plotted in Fig. 16. 8[ 
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Figure 6.8: Invariant mass of eight gamma quanta from the rj meson decay. Sohd 
hne correspond to Gaussian function fitted to the histogram. 

By fitting the Gaussian function to the histogram, a resolution of rj meson invariant 
mass reconstruction in pp — )■ ppr] — )■ ppin^ — >■ ppS'y reaction with WASA-at-COSY 
detector is determined to o" = 31 MeV/c^. 

Present software setup and restrictions imposed during analysis allow to recon- 
struct r] meson only for 3.1% of all pp ppr] ppA-n^ pp8j events. 



7 Time of measurement 



The main goal of this thesis is to estimate time of measurement, with the 
WASA-at-COSY detector, for the reaction pp — )■ pprj — )■ pp47r° — )■ pp87 allowing 
to improve the present knowledge about the BR{ri 47r°). Up to now only an 
upper limit of the branching ratio has been known ^4] : 

BR{7] 47r°) < 6.9 • 10"^ CL = 90%. (7.1) 

which is determined by Crystal Ball Collaboration [5]. 

In general a number of reconstructed 1] — 47r° decays can be written as 

N^^AnO = BR{r] 47r°) ■ ■ ■ L ■ (7.2) 

where N^^^^^o and BR stands for number of observed events and branching ratio 
respectively, is acceptance of the detector, L denotes luminosity, cr^ corresponds 
to cross section for rj meson production and At is time of measurement. 
Rewriting it as an equation for branching ratio, we find 

Statistical errors of branching ratio depends on uncertainty of all variables 
(cr(A^^_j.4^o), a{Aj^), o"(cr^), (t(L) and o"(At)). However, uncertainty related to A^, a^, 
L and At can be assumed to be negligibly small with comparison to c'"(A^^_^4^o). 
Under this assumption a{BR) may be expressed as: 

can be established using Monte Carlo simulation as a ratio of observed events 
to produced events, = N^^ detected/ produced- In the determination for a number 
of detected a missiug mass technique is applied. The missing mass technique is 
essential to distinguish the decay of rj meson from prompt production of 47?*^ in 
proton-proton collision (main background of — )■ 47r°). In Fig. 17.11 distribution of 
missing mass reconstructed from sample of simulated pp — )■ ppr] — )■ ppAn^ (solid, 
black line) and its main background pp — )■ ppAn^ (dashed, red line) is plotted. These 
events obey all restrictions mentioned above in section[6l Full width at half maximum 
(FWHM) for 7] signal is approximately 3 MeV. Taking as A^,j^4^o a number of entries 
in the range of 6 MeV (iFWHM) around the mean value one obtains that is 
equal to about 1.5%. 
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Figure 7.1: Missing mass for the pp — ppX reaction. Solid, black line denotes re- 
action where the rj meson is produced. Dashed, red line denotes prompt production 
of 47r° in proton-proton collision. Both histograms are results of Monte Carlo sim- 
ulations of 10^ events for pp pprj — )■ ppAir'^ — )■ ppS'j and pp — i- ppAn^ — )■ pp8'~f 
reactions with Tbeam = 1-3 GeV. 

In order to determine uncertainty of number of observed events, information 
about background reactions is needed. The value of N^^^^o can be in principle 
calculated as: 



■^Tj^i-K'^ -^observed -^background \-^r]^4iT^^ ~l~ ^background) ^background' C^"^) 

Under assumption that events observed in experiment can be described by Poisson 
distribution and that Nobserved as well as Nbackground can be measured independently 
and that N^j^/^^o is negligible compared to Nbackground, uncertainty of N^^/^^^o equals 
to: 

(^{Nri^Ano) = a/2 ■ Nbackground (7.6) 

Using equation 17.21 to describe Nbackground we find: 

(j(A^^_^47ro) = a/2 ■ A^j^o ■ a^T^o ■ L ■ At, (7.7) 

where L and At denote the same variables as in equation 17. 2[ A^^^o stands for ac- 
ceptance of the detector for prompt four pions production and a^-^o corresponds to 
cross section for pp — > ppAir'^ reaction. Employing equation 17.61 in equation 17.41 one 
obtains: 
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^47r0 


■ a4T,o 




■ . 


L-At 



(7.8) 



Using the same procedure for ^47^0 as it is described above for A^, one can obtain 

^4^0 = 0.5% 

Another variable in equation 17.21 which has not been mentioned yet is cross 
section. The cross section of the 7] production (cr^) is taken from the parameterization 
of experimental data dB [HI [131 [H [15] • 

To describe experimental cross section for the 1] meson production in proton- 
proton collision, third order polynomial {a + bQ + cQ^ + dQ^) is fitted to data. 
The result of fitting is shown in Fig. 17.21 The values of polynomial coefficients are 
presented in Tab. 17.11 




Figure 7.2: Cross sections for the 77 meson production in proton-proton collision. 
Data are from Refs. [TT], [T21 [13, [H [U] and the superimposed line indicates result 
of the fitted polynomial a + hQ + cQ^ + dQ^. 

The cross section for the prompt pions production is taken from parameterization 
14.41 (see section Hj) and its shape is shown as a green line in Fig. 14.21 
Time {At) is a parameter. 

Using equation 17. 8[ statistical error of branching ratio as a function of time can be 
plotted and it is shown in Fig. 17.31 under assumption that luminosity is constant and 
equal to L = 10^^ cm~^s~^. The result is obtained for proton-proton collision with 
beam kinetic energy of Tbeam = 1-3 GeV . Upper limit of a branching ratio for the 
present value of the decay 77 — )■ 47r° (17.11) is marked and it can be achieved after 13.3 
hours of measurement using WASA-at-COSY detector. However, this result depends 
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Coefficient 


Value 


a 


-0.0164 ± 0.0040 


h 


177.3 ±4.4 


c 


-55 ± 17 


d 


-2.3 ±6.2 



Table 7.1: Coefficients of polynomial {a + bQ + cQ^ + dQ^) fitted to cross sections 
for the 7] meson production in proton-proton collision as shown in Fig. 17. 2[ Q is 
expressed in units of GeV. 
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Figure 7.3: Statistical error of branching ratio as a function of time. Blue square 
presents value of upper limit of branching ratio ^ . Black line shows value of statis- 
tical error which can be achieved using WASA-at-COSY detector for measurement 



at Tf 



beam 



1.3 GeV with L = 10^^ cm-^s"^ and time At. 



strongly on the beam energy. Dependence on the excess energy is hidden in cross 
sections: 0-4^0 = cr^noiQ), o-^ = (JniQ) and acceptances: A^j = A^^^Q), A^^o = A^^o{Q). 
Thus, by rewriting equation 17.81 for time of measurement as a function of Q, we find 

A^iQ) ■ L ■ a^BK) ■ a^^iQY ^ ' ^ 

In order to find the optimum energy for the studies of the pp — )■ ppr] — )■ ppAn^ 
reaction, a time of measurement is plotted as a function of excess energy (Fig. 17.41) 
with fixed value of a{BR) = 6.9 ■ 10"^. 

Fig. 17.41 indicates that time of measurement needed to achieve a given accuracy 
of BRij] — )■ 471°) decreases with smaller values of excess energy for the t] meson 
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Figure 7.4: Relation between time of measurement needed to achieve statistical 
uncertainty of 6.9 • 10~^ and excess energy for the rj meson production in proton- 
proton coUision assuming L — 10^^ cm~'^s~^. 



production in proton-proton collisions. 
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8 Summary and conclusions 



The main aim of this thesis was to estimate the time of measurement of 
1] — > 47r*^ decay using WASA-at-COSY detector for which the present branching 
ratio (BR) limit can be improved. To this end a Monte Carlo simulations of 
PP ~^ PPV ~^ PP^'k'^ ppS'y reaction kinematics and simulations of WASA-at-COSY 
detector response were performed. 

The first step of signals analysis was the estimation of how many of rj decays 
could be reconstructed. For this purpose geometry and reconstruction efficiency of 
WASA-at-COSY detector for registration of gamma quanta from rj — > 47r° — > 87 
decay were examined. The obtained result shows that only about 3% of 77 — > 47r° — > 
87 events can be reconstructed. Resolution for reconstruction of eight gamma quanta 
invariant mass was established to 31 MeV/c^ (standard deviation). 

In the next step of analysis the detector acceptance for the pp — )■ ppr] — )■ ppAir^ — )■ 
ppS'j reaction {A^ = N,i deiected/N^ produced) was estimated to = 1.5%. A number 
of A^^ detected Gvcuts was extracted from the missing mass distribution of pp ppX 
reaction. 

The last step of analysis was to evaluate time of measurement to obtain statistical 
uncertainty equal to the current upper limit of BR. For proton beam with kinetic 
energy of Tbeam = 1-3 GeV and luminosity of L = 10^^ cm~'^s~^ the obtained time 
is equal to about 13 hours. In addition, based on the parameterization of cross 
section for pp — > pprj and pp — > ppAir^ reactions the time needed for measurement 
was also studied as a function of excess energy for 77 meson production in proton- 
proton collision. From the determined dependence it was concluded that the best 
conditions for the registration of the pp — >■ ppr) pp47r° — )■ ppS'-f reaction with the 
WASA-at-COSY detector are for beam energies close to the kinematical threshold 
for pp — )■ pprj reaction. 

Obtained results indicate that statistical uncertainty at level of present branching 
ratio limit is achievable with the WASA-at-COSY experiment. However, its signif- 
icant improvement will require long term measurement. In order to improve the 
current limit by about one order of magnitude about 41 days of the measurement 
at Tfeeam = 1-3 GeV and luminosity of L = 10^^ cm''^s~^ is needed corresponding to 
the production of about 10^ rj mesons. 

Up to now almost 4-10^ decays of the rj meson have been collected using 
WASA-at-COSY detector. 

Consequently, the outcome of this thesis demonstrates that it is worth to search 
in these data for the rj — >■ 47r° decay. 
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A Kinematical Fit 



Distribution of rj meson mass reconstructed as invariant mass of eight gamma 
quanta is quite broad (see Fig. 16.81) . In order to improve this resolution, kinemati- 
cal fit technique is applied. Kinematical fit varies measured quantities for particles 
within the measurement uncertainty to match constraints of working hypothesis and 
optimize precision of the studied variables. In this case, as constraints it is imposed 
that invariant mass of every pair of gamma quanta matched by routine (see sub- 
section 16. 2p should be equal to neutral pion mass. Quality of fit is measured by 
calculation. Varied quantities for protons and gamma quanta are kinetic energy E^^, 
polar angle 6 and azimuthal angle 0. Uncertainty of these quantities are obtained 
on a basis of the Monte Carlo simulations by the comparison of distribution of the 
difference between true and reconstructed values {qmc ~ Irec)- As an example of this 
kind of distribution, kinetic energy of gamma quanta is shown in Fig. lA.li 
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Figure A.l: Difference between true and reconstructed kinetic energy of gamma 
quanta determined from signals in the electromagnetic calorimeter. In order to de- 
termine standard deviation, Gaussian function is fitted. 

All standard deviations obtained from Gaussian function fitting for each variable 
are presented in Tab. lA.li 

Kinematical fit for input data mentioned above did not work well. Therefore, 
efficiency of gamma quanta matching algorithm is performed. For the best set of 
gamma quanta according to the routine, a distribution of number of properly as- 
signed gamma quanta is presented in upper histogram in Fig. IA.2[ As one can notice 
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Appendix Kinematical Fit 
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Table A.l: Standard deviations used in kinematical fit program, 
only in about 1/3 of reactions gamma quanta are fully properly assigned. 
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Figure A. 2: Number of gamma quanta correctly matched into pairs by the used 
algorithm. The upper histogram shows gamma quanta for the best set chosen by 
the algorithm. The bottom one presents the choice with the biggest number of 
properly assigned gamma quanta from four best possibilities found by the matching 
algorithm (see subsection 16.21) . 

For comparison purposes, beside the best set of gamma quanta according to 
the algorithm, three additional sets are checked. Additional sets are these with 
a larger sum of differences between pion mass and masses of matched pairs (see 
subsection 16.21) . The chosen set taken into account in Fig. IA.2l is the one for which 
number of correctly matched gamma quanta is the biggest. The result of this test is 
shown in bottom histogram in Fig. IA.2I One can observe a significant improvement 
in gamma quanta assignment. 
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Signals reconstructed from Monte Carlo simulation can be compared with true 
values, however in experiment other criterion of properly matched gamma quanta is 
needed. In experiment value from the kinematical fit can be used. 

Similar comparison to this above is done for events with value as criterion 
of proper gamma quanta assignment. In Fig. IA.3| right hand side histogram shows 
invariant mass of eight gamma quanta after kinematical fit taking into account only 
the best outcome of the matching algorithm. Histogram in the left panel shows result 
after kinematical fit taking the outcome with the best out of four best possibilities 
found by the matching algorithm. Set assumed to be the best is the one with 
value closest to the expected value of distribution of pp — )■ pprj — )■ ppAn^ — )■ ppS'y 
reaction. 




Figure A. 3: The right hand side histogram shows invariant mass of gamma quanta 
for the best set of the gamma quanta found by the matching algorithm. Histogram 
on the left hand side presents the choice with the value closest to the expected 
value of distribution of pp — )■ ppi] — )■ ppAn^ — )■ ppS'j reaction. 

FWHM of the distribution in Fig. |X3] (left) is equal to 40 MeV/c^ (a = 
17 MeV/c^) and it is almost two times smaller compared to one obtained in sub- 
section 16.41 (cr = 31 MeV/c^). One can see that number of entries is more than 
two times higher in case when four sets of gamma quanta are checked. However, 
still low statistics and non-gaussian distributions of both invariant mass histograms 
in Fig. IA.3I indicate that kinematical fit does not work satisfactorily for that input 
data. Therefore, further studies will be continued in the future. 
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